A simple model for self organization of bipartite networks 
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We suggest a minimalistic model for directed networks and suggest an application to injection 
and merging of magnetic field lines. We obtain a network of connected donor and acceptor vertices 
with degree distribution and with dynamical reconnection events of size As occurring with 

frequency that scale as 1/As^. This suggest that the model is in the same universality class as the 
model for self organization in the solar atmosphere suggested by Hughes et al. 0. 

PACS numbers: 05.40.-a, 05.65.+b, 89.75.-k, 96.60.-j, 98.70.Vc 



In a number of physical systems one observes emer- 
gence of large-scale structures, caused by growth of small- 
scale fluctuations. For example, 1) the energy flows from 
small to large scales in 2-d turbulence, 2) the matter dis- 
tribution in the universe is highly inhomogeneous in spite 
of a presumably uniform energy distribution at its origin, 
and 3) the magnetic field lines reconnection and sunspot 
activity is able to generate solar flare activity with burst 
sizes that by far exceed excitations associated to the indi- 
vidual convection cell on the solar surface. In fact, often 
the emerging large-scale structures exhibit scale-free fea- 
tures over substantial range of scales, as e.g. the sun 
spots 0, and solar flare activities H, ^ |1| . 

Recently it has been realized that many complex 
networks exhibit scale-free topologies 0, , including 
in particular the topology of sun spots connected by 
magnetic field lines [ij, |3] ■ In general, the first theoretical 
framework for emergence of power law distributions 
was the Simon model j^], featuring a "rich get richer" 
process, that recently has been developed into preferen- 
tial attachment to explain scale-free networks An 
alternative approach to generate large-scale features 
from small-scale excitations is provided by the self 
organized critical (SOC) models 0, Ull UJ which m 
their traditional versions propose a scenario for the 
fractal pattern of activity that is observed in systems 
with extreme separation of timescales. Hughes et al. 
has proposed a SOC like mechanism for cascades of 
reconnection of magnetic field lines in the solar atmo- 
sphere, using a plausible number of processes associated 
to diffusion of sun spots and reconnection of crossing 
fleld lines. In this paper we suggest a simpler model, 
assuming only two processes, merging and creation, in 
an on going dynamics of vertices connected in a network. 



We flrst review the basic process of merging-and- 
creation ( originally proposed by 0, Q) in a formu- 
lation that is closest to the network interpretation which 
we will discuss later. The model describes the evolution 
of a system of many elements i = 1,2, N that each is 
characterized by a scalar qi that may be either positive 



or negative. One may think of the scalar as a helicity or 
as a quantiflcation to which extent an element/ vertex is 
a donor or an acceptor. The model describes a situation 
in which the elements in the system redistribute their 
respective charges qi according to 
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With {k,l) selected independently from {i,j) these two 
processes define one of the many possible realizations of 
the model. Other realizations include different combina- 
tions of correlations between {k,l) and {i,j). For exam- 
ple, one may select k — j and I — i. For any choice the 
obtained scaling is as reported in Fig. ^ 

The main features obtained numerically are presented 
in Fig. 2] Fig. ^a) illustrates the steady state after a 
transient time ~ N updates per element, starting from an 
initial "vacuum" with qi = 0,Vi = l,2,...iV. The figure 
shows the extreme range of g-values at any time. The 
subsequent dynamics of the extremes is also reflected in 
the trajectory of a winner-element which, when merged, 
is re-identified as the merged element. One observes that 
this winning element exhibits an intermittent dynamics 
with size-changes Aq of all magnitudes. The distribution 
of these changes as well as a wide set of other properties 
is in fact scale invariant. The cumulative distribution of 
g-values, Fig. ^b), is a scale-free distribution. 
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with 7 = 2. With asymmetric initial condition, say 
qi = 10, Vi = 1,2,...A^, as illustrated in Fig. ^c), the 
system self-organizes by concentrating all of the initial 
asymmetry to one of the elements. All other elements 
are distributed in exactly the same way as with the "vac- 
uum" initial condition (compare Figs.^b) and^Jc)). 

Fig.njd) shows the distribution of changes in Ag under 
steady state conditions. There are two possible ways to 
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1: Main features of the basic model. In panel a) we 
show the development of the elements with largest positive 
Qmax and largest negative qmin together with a particular tra- 
jectory, where we always follow the "winner" in the merging 
process (solid line). The simulation is done for a system of 
size A'^ = 10* elements, and the time-count is in updates per 
element, b-e) Cumulative plots of steady state properties, b) 
The size distributions of the positive and negative q for 3 dif- 
ferent system sizes, c) The size distribution when we start 
a system of size N = 10^ with initial condition qi = 10, Vz. 
One observes that all excess q is moved to a single element. 

d) Two variants of a histogram of size changes. Aqmm is 
defined by following the winner, meaning that we plot the 
difference in size between the largest of qi,qj before merging 
and the merged unit after the merging. Ag^m corresponds 
to the change from any of the two qi, qj to the merged unit. 

e) The size of changes defined as the losses of absolute q in 
merging events where qi and qj are of different signs. 



characterize these changes. One may quantify them by 
considering the difference between the merged element 
ili + Qj) f-iid s-iiy of the two qi or qj merging elements. 
In that case one observes a cumulative distribution for 
changes P{> Aq) = l/(Aq), (compare full drawn Hue 
and dashed line in Fig. ^d)). This distribution closely 
resembles the overall distribution of q values. Alterna- 
tively one may quantify the dynamics by following the 
winner at each merging, thus defining the Aq as the dif- 
ference between the largest q before and the largest q 
after the merging. In that case one expects the probabil- 
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which with 7 = 2 from Fig.^b) predicts exponent 3 veri- 
fied by simulations, see Fig.^d). Finally Fig.^e) shows 
the size-distribution of annihilation events, defined as 
events where two elements of different signs merge. The 
distributions of these annihilations are governed by the 
same considerations as in Eq. 01 and accordingly scales 
with exponent t — 3. 

Now we explore the reason for the 7 = 2 scaling be- 
havior. We consider the version of the model where one 
excites the system by randomly picking a zero element 
and assigning it a +/~ value: 



(5) 
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were r is a random number picked from a symmetric nar- 
row distribution F{q). This update is one of many pos- 
sible versions that all produce the same scaling results 
as shown in Fig. 1, we here consider it because it is the 
simples to treat analytically. The differential equation, 
describing the evolution of the model reads llJ] 



dPjq) 
dt 



// dq2dqiS{q - qi - q2)Piqi)Piq2) 

J J-OG 

2P{q)+F{q), (7) 



which have been shown to give a steady state distribu- 
tion with the asymptotic behaviour P(g) oal/q^ For 
pedagogical reasons we here present an alternative solu- 
tion, that also opens for some insight into the amazing 
robustness of this model. In terms of the Fourier trans- 
form p{ijo) — J dqe^^'^'^ P{q) the steady state equation is 



(8) 



The important property is that p{uj) — 1 oc — |w| for small 
uj. A positive creation probability F(q) with a finite sec- 
ond moment ensures this which leads to P{q) oc q^^ for 
large q. Thus the exponent 2 will be a common property 
for a large class of variations of the basic merging and 
creation mechanism. As an example F{q) — exp(— |g|)/2 
gives 



P(q) = ---Re[5n(*g)] = ^- A 



(9) 



where Sii{z) is a Lommel function and P{k = 0) = I/tt. 
Also the localization of positive excess {q)N can be un- 
derstood, since a symmetric F{q) implies an even contin- 
uum solution and thus that all excess will occupy a zero 



3 




10 
-10 



1 ..... 1 1 i. . . . 






■ ir 1 





25 



50 



75 



FIG. 2: Merging and creation in 1-dimension: At each time- 
step one selects a coordinate i between 1 and N=100. If qi = 
0, one create a +/— pair at position i and one of its neighbors. 
If 5i 7^ then one moves qi one step to either left or right and 
adds it to the q aheady present at that position. In upper 
panel one sees the time evolution over in total 100 updates 
per site, whereas the lower panel is a snapshot of the final 
configuration. One may notice that merging of light (positive) 
and dark (negative) sites leads to annihilation of both. The 
steady state distribution of q is scale invariant with the same 
exponent 2 as in the basic (infinite dimensional) model. 



q measure around go ~ {q)^- For ^ discrete simulation 
this means a single q as illustrated in Fig. QJ;. 

Finally, for application in real physical situations, 
it is also of interest to explore the behaviour of the 
merging-and-creation scenario in finite dimensions. As 
was reported by 0, Q, 0| then the observed scaling 
is robust, even when we confine the elements to 
diffusive motion in 1 dimension, provided that creation 
of +/- pairs occur close to each other. For a visualization 
of the dynamic behaviour we in Fig. [21 show the evolving 
system in 1-d. 



We now consider a network implementation where each 
element is a vertex and its sign corresponds to the number 
of in- or out-edges. Thus the above scenario is translated 
to a network model in which donor (q > 0) and acceptor 
{q < 0) vertices are connected by directed edges, see Fig. 
|3Ja) and (b). Each vertex may have different number 
of edges, but at any time a given vertex cannot be both 
donor and acceptor. Further, in the direct generalization 
of the model, we allow several parallel edges between any 
pair of vertices. At each time-step two vertices i and j 
are chosen randomly. The update is then: 

• Merge the two random vertices i and j. There are 
now two possibilities: 

a) If they have the same sign all the edges from i 
and j are assigned to the merged vertex. Thereby 
the merged vertex has the same neighbors as i and 
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FIG. 3: The network realization of the model, a) and b) 
illustrate possible merging moves. Positive vertices (donors) 
are vertices with outgoing edges and negative (acceptors) with 
incoming edges, c) The dynamics of the average number of 
edges per node, (E), (upper curve) and the average number 
of neighbors {nn) (lower curve), A*' — 10*. d) The cumu- 
lative probability distributions, N — 10^, for: number of 
edges incoming or outgoing from a node, E (solid curve); 
number of neighbors, nn (dotted curve); edge density, ps 
defined as the number of parallel edges connecting two ver- 
tices (dashed curve). The distributions for all quantities are 
scale-free P{>s) ~ l/s''"^ with 7 = 2. e) The cumulative 
probability distributions for the changes in number of edges 
due to merging, AE and number of neighbors Ann. The dis- 
tributions are power-law P{> As) ~ As^~'^ with exponent 
r = 27 - 1 = 3 from Eq. 4. 



j had together prior the merging, see Fig.|2fa). 
b) If i and j have different signs, the resulting ver- 
tex is assigned the sign of the sum qi + qj . Thereby 
a number max{\qi\, \qj\} — \qi + qj\ of edges are an- 
nihilated in such a way that only the two merging 
vertices change their number of edges. This is done 
by reconnecting donor vertices of incoming edges to 
acceptor vertices of outgoing edges, see Fig. (Sfb). 

• One new vertex is created of random sign, with one 
edge being connected to a randomly chosen vertex. 

On the vertex level this network model can be mapped 
to the above model for merging and creation, and thus 
predict similar distributions of vertex sizes, as seen by 
comparing solid line in Fig.l^Jd) with Fig.^b) and distri- 
butions of annihilations in Fig.|3Ie) and Fig.nje). How- 
ever, the network formulation provides additional insight 
into the excitation process that drives the whole distribu- 
tion. That is, starting with a number of empty vertices 
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Qi = 0, the creation process generates vertex antivertex 
pairs on small scale which subsequently may grow and 
shrink due to merging and creation as illustrated in Fig. 
|3Jc). One can see, that when the system has reached the 
stationary state, the average number of neighbors (nn) 
is nearly constant with small fluctuations while the fluc- 
tuations in the average number of edges, (i?), are much 
larger. Further one notices that the evolution of {E) is 
asymmetric, in the sense that increases are gradual, while 
decreases are intermittent with occasional large drops in 
(E) . These drops primarily correspond to the merging of 
vertices of different signs, where a large number of edges 
may be annihilated. This process is quantified in Fig. 
He). 

The network model opens for a new range of power- 
laws [Bj associated to the connection pattern and dynam- 
ics of reconnections between the vertices. In this connec- 
tion it is interesting that the number of edges per vertex, 
E, is distributed with scaling P{E) oc 1/E'^. This was 
also obtained for the "number of loops at foot-point" in 
[3- In addition, the distribution of reconnection events 
Pa(A£') oc 1/AE^ is distributed as the "flare energies" 
in the model of Ref. In our model the event size 

is simply the change in the number of edges (AE) when 
two vertices merge which gives the exponent -3 as shown 
in Eq. [3 In the model of Hughes et al. .5] the event size 
is a more complex quantity related to cascades of cross- 
ings of field lines, and the energy release is associated to 
the number of lines that thereby decrease their length. 
The non trivial fact that we obtain the same exponent 
suggests that the two models are in the same universality 
class, which means that our minimalistic model captures 
the main features of a presumably much larger class of 
more detailed and realistic models. 

Also we would like to mention that the distribution of 
the number of parallel edges for connected pairs of nodes 
is also scale invariant P{>pe) oc p]^'^ see Fig.|SId). This 
illustrates robustness of the mechanism: The dynamics 
of merging vertices appears very different when it is 
viewed from the "dual" space of tubes of edges between 
vertices, p^, nevertheless the same exponent 7 = 2 is 
obtained. 



In conclusion we have discussed a new mechanism for 
obtaining scale-free networks of connected donor and 
acceptor vertices. The model predicts power-laws of 



node degrees with a 1/s^ distribution, and of recon- 
nection events with a 1/As'^ distribution. The scenario 
thus provides a generic framework to generate networks 
with large-scale features from small-scale excitations un- 
der steady state conditions, and may thus complement 
preferential growth which provides scaling only under 
persistently growing conditions 0. Viewed as SOC, 
the merging-creation scenario provides "scaling for free" 
in the sense that it is robust to multiple simultaneous 
updates. The key process of both constructive (equal 
sign) merging and destructive (opposite sign) merging 
[l7| should be an important ingredient in a number of 
dynamic systems, and in particular appears to be ap- 
pealing minimalistic model with possible connection to 
reconnection and creation of solar flares. 
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